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We discuss uncertainties and possible sources of errors associated with the determination of the diffuse Galactic 7-ray 
emission using the EGRET data. Most of the issues will be relevant also in the GLAST era. The focus here is on issues that 
impact evaluation of dark matter annihilation signals against the diffuse 7-ray emission of the Milky Way. 



1. INTRODUCTION 

Diffuse emission (DE) from the Milky Way dom- 
inates the 7-ray sky. About 80% of the high-energy 
luminosity of the Milky Way comes from processes 
in the interstellar medium (ISM). The DE traces in- 
teractions of energetic particles, primarily protons 
and electrons, in the ISM, thus delivering informa- 
tion about cosmic -ray (CR) spectra in distant loca- 
tions LLJ. The DE may contain signatures of exotic 
physics, e.g., interactions of dark matter (see (D for 
references). Calculation of the DE requires first cal- 
culating the CR spectra throughout the entire Galaxy 
llJ). The components of the DE model (7r°-decay, in- 
verse Compton - IC, bremsstrahlung) are not inde- 
pendent and can not be arbitrarily re-scaled (e.g., as in 
lis]). The DE is the celestial foreground for the study 
of 7-ray point sources and the extragalactic DE which 
may contain information about the early universe. 

A comparison between the EGRET data and the 
Galactic DE model used by the EGRET team reveals 
an "excess" above 1 GeV of approximately a factor 
of 2 in all sky directions 1214 151 . Therefore, we need 
to look at the uncertainties associated with both the 
data and the model to understand the origin of this 
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excess. In this paper, we consider the sources of sys- 
tematic uncertainties in the EGRET calibration, data 
handling, and in models of the DE. 

2. SOURCES OF UNCERTAINTIES 

2.1. EGRET data analysis 

EGRET was in orbit between April 1991 and June 
2000, and operated for much longer than its 2-year de- 
sign lifetime. Owing to ageing of the spark chamber 
gas and the limited number of gas changes possible 
during the mission, the best data taking occurred dur- 
ing the early years of operation. Details of pre- and 
in-flight calibrations are given in [6 7 1. 

The data that we discuss here were taken during 
Cycles 1-4, the "prime" years of EGRET operation. 
These are the data that were used to produce the 3rd 
EGRET catalog of 7-ray point sources (8), and so for 
these data subtracting the (generally variable) point 
source contribution to the DE can be done accurately 
and correctly using the catalog fluxes. The overall ex- 
posure, accumulated during the individually pointed 
observations, is rather uneven |7|. 

The DE model used for the source detections 
was derived on the assumption that the CR proton 
and electron spectra do not change shape across the 
Galaxy, and that the CR density is proportional to the 
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surface density of gas . This leads to incorrect pre- 
dictions at intermediate and high latitudes. 

The calibration of an instrument in orbit is a com- 
plicated matter. Issues that cannot be controlled may 
affect the instrument performance, for example the 
launch vibrations and accelerations, while radiation 
exposure and unknown factors occur during the whole 
mission lifetime. Several failures occurred during 
egret's operational years: the spark chamber read- 
out (switched to the redundant side), the gas circula- 
tion pumps, and some PMTs of the trigger telescopes. 
This complicates the analysis since the symmetry of 
the instrumental response was lost, and thus the direct 
comparison to the pre-flight calibration. 

The efficiency of EGRET for registering recon- 
structable events from 7-ray detections decreased 
with time as the spark chamber gas aged. Several 
times during the mission the gas in the spark chamber 
was replenished to improve the degraded efficiency. 
Efficiencies were determined based on normalization 
to "standard candles" (pulsars) where possible, or by 
comparing to the DE predictions |7|. The individual 
EGRET viewing periods were corrected for the vary- 
ing spark chamber performance under the assumption 
that the energy dependence of the correction did not 
vary with time |7|. This results in time-dependent 
systematic uncertainties. 

Several other uncertainties also must be considered 
fTl : contribution of unresolved point sources; fish- 
eye effect corrections for events under large viewing 
angles; contamination by Earth albedo events |9|; a 
direction-dependent PSF which results in residuals in 
source-subtracted DE maps; several orbital reboosts 
which changed the particle background of the detec- 
tors; plus, late in the mission, operation in different 
modes (narrow, strip) for which the instrument per- 
formance is less well understood. 

After these considerations, any model of DE must 
be convolved with the instrumental PSF before com- 
parison with the data. The effect is large below 1 GeV 
affecting the overall spectral shape. This is especially 
important for small areas on the sky. For comparisons 
with data obtained by a different instrument, a de- 
convolution of the data is preferred |2|: in this case 
the procedure is model-dependent since the predicted 
spectrum is convolved with the PSF first, and subse- 
quently the appropriate factors are applied to the data. 

The EGRET team made corrections for the de- 



scribed effects taking into account the limited statis- 
tics. The quoted error bars |7 j are not separable into 
statistical and systematic parts (e.g., as in JSJ) because 
of the implicit time dependence and interplay of the 
applied corrections. 

2.2. CR propagation 

The Galactic diffuse 7-ray emission is the product 
of particle interactions with gas in the ISM and the 
interstellar radiation field (ISRF). Therefore, its spec- 
trum depends on the particle spectra throughout the 
entire Galaxy (T]. Diffusion, energy losses, and other 
processes change the particle spectra during propa- 
gation in the Galaxy over long time periods. Addi- 
tionally, the CR isotopic composition changes. Dur- 
ing propagation within the heliosphere, solar modu- 
lation further alters the particle spectra. These spec- 
tra are what can be measured directly using satelUte- 
and balloon-borne instruments. To calculate the 7-ray 
emission we thus need to "reconstruct" the particle 
spectra in the ISM. 

Gamma rays from 7r"-decay are produced in the 
same processes as other secondaries such as antipro- 
tons and positrons, which can also be used to test the 
models of the DE [2' 101. A realistic model of CR 
propagation has to include the gas and source dis- 
tributions, ISRF, nuclear and particle cross sections 
and nuclear reaction network, 7-ray production pro- 
cesses, energy losses, and solve the transport equa- 
tions for the individual CR species. The propaga- 
tion parameters, the diffusion coefficient, halo size, 
Alfven speed, convection velocity, and so on, which 
depend on the assumed propagation model, are de- 
rived by comparing model predictions with CR data, 
such as secondary /primary nuclei ratios, e.g. B/C, and 
radioactive isotope ratios like ^^Be/^Be. Therefore, 
the accuracy of the nuclear cross sections is one of 
the major concerns for uncertainties in all subsequent 
analysis tasks. 

Semi-empirical estimates of the nuclear cross sec- 
tions 1 11 12 1 have been tuned to match the isotopic 
production cross sections between '^400 and 700 
MeV/n. They can not provide the same accuracy 
over the whole energy range which is required for CR 
propagation calculations. Furthermore, the cross sec- 
tions often have resonances at energies below '^200 
MeV/n. These are particularly important for propaga- 
tion models including diffusive reacceleration. This 
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adds to the overall uncertainty in determination of the 
propagation parameters (e.g., see |13|). 

The interpretation of the peak in the B/C ratio de- 
pends on the assumed propagation model. Different 
models predict different rigidity dependencies for the 
diffusion coefficient. Therefore, determination of the 
diffusion parameters is model-dependent. Two exam- 
ples are the reacceleration |14| and convection |15| 
models. Reacceleration is a Fermi 2nd order accel- 
eration in the ISM; it predicts a diffusion coefficient 
rigidity power-law index of 1/3. Convection is essen- 
tially a Galactic wind moving away from the Galactic 
plane which affects mostly particles below a few Ge V; 
consistency with high-energy data requires a diffusion 
coefficient rigidity index of 0.6. A recently developed 
"damping" model includes CR particle - MHD wave 
interactions which lead to a concave shape in the dif- 
fusion coefficient, where the true rigidity dependence 
is derived via self-consistent modeling 1 16 1 . 

Above 1 GeV, Tr'^-decay and IC scattering are 
the dominant 7-ray production processes. The tt"- 
production cross section is based on data from the 
1960s that themselves have large error bars. The pro- 
posed parameterizations ( 1171 and references therein) 
fit these data well, but they do not provide the required 
accuracy. New details are still being added, e.g. the 
diffraction dissociation and violation of the Feynman 
scaling 1 18 1, while various Monte Carlo event gener- 
ators produce results which differ by ^20-30% 1 19 1. 
The IC scattering contribution is dependent on the un- 
derlying model of the ISRF, and the treatment of the 
scattering process using the proper ISRF angular dis- 
tribution |20|. Conventionally, the ISRF is assumed 
to be isotropic, which is only true for the cosmic mi- 
crowave background. Using the full angular distribu- 
tion for the ISRF will alter the model predictions 1 20 1 . 

The random nature of CR sources leads to fluctua- 
tions of CR intensity in space and time. Very high en- 
ergy (VHE) electrons cannot propagate far from their 
sources due to rapid energy losses. Their interstel- 
lar spectrum is expected to exhibit imprints from their 
nearby sources 1211 . 

Finally, "solar modulation" changes the interstellar 
spectra of CR particles below ^20 Ge V/n during their 
propagation within the heliosphere. It is a combina- 
tion of the effects of convection by the solar wind, dif- 
fusion, adiabatic cooling, drifts, and diffusive acceler- 
ation. The theory of solar modulation is far from com- 



plete [22]; current models are based on the solution of 
Parker's transport equation. The unknown is the in- 
terplay of the different terms and the CR spectrum in 
the ISM. Spherically symmetric solutions, the force- 
field 1 23 1, and Fisk |24| approximations are most of- 
ten used. However, they include only the effect of 
adiabatic losses. The Pioneer, Voyager-', and Ulysses 
missions contributed to understanding the global as- 
pects of modulation. 

2.3. Astrophysical input 

Most of the interstellar gas is in the form of neu- 
tral hydrogen H I and H2. The H2 gas is not observed 
directly but via the mm-wave spectral line of a rota- 
tional transition of CO. A proportionality is assumed 
between CO surface brightness and H2 column den- 
sity. H I gas is detected by its 21 cm line. Deter- 
minations of H I column densities rely on difficult- 
to-verify assumptions about the spin temperature of 
the gas, and in many directions colder H I gas is seen 
in absorption against warmer background gas, mak- 
ing overall distributions of H I difficult to disentangle. 
The spatial distribution of the gas is derived using ve- 
locity measurements and assuming a rotation curve 
for the Milky Way. However, streaming motions and 
velocity dispersions between individual gas clouds in- 
troduce systematic uncertainties. In the outer Galaxy 
the gradient of velocity approaches zero at large dis- 
tances. In the inner Galaxy the "kinematic distance" 
is double valued except toward the Galactic center 
and anticenter, where the line-of-sight velocities pro- 
vide no useful kinematic information. Therefore, no 
unique solution exists for the 3D distribution of gas. 
Studies starting with the same data produced gas dis- 
tributions which are different in many details 14 1261 . 

The CR source distribution is also not very well 
known. The pulsar distribution |27J is too peaked 
to reproduce the CR gradient in the Galaxy. Con- 
sistency with the latter requires the X-factors (essen- 
tially H2/CO ratio) to increase with the distance from 
the Galactic center II28I . Such behavior indeed is ob- 
served in other galaxies | 291. 

The ISRF modeling incorporates a stellar distribu- 
tion model, a model for the dust distribution and prop- 
erties, and a treatment of scattering, absorption, and 
subsequent re-emission of the stellar light by the dust 

•^The Voyager 1 spacecraft is currently approaching the outer 
boundary of the solar system 1251 . 
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1301 . The stellar emission model incorporates details 
for the distribution of stellar types in discrete geo- 
metrical components (bulge, disc, spiral arms, etc.), 
adjusted to agree with results from experiments such 
as 2MASS, SDSS, and others. Dust is modeled with 
a mixture of PAH, graphite, and silicate. Details of 
the dust absorption and scattering efficiencies, abun- 
dances, and size distribution are included in the scat- 
tering and heating calculations. The dust is assumed 
to follow the Galactic gas distribution and metallicity 
gradient. Uncertainties in each of these inputs con- 
tribute to the overall uncertainty in the ISRF model- 
ing. The model can be compared with the ISRF only 
at the solar system position. The deviations from an 
earher model 1513 II give an idea of possible system- 
atic errors. 

3. DISCUSSION AND CONCLUSION 

We discussed only the most obvious sources of un- 
certainties; the systematic effects are numerous, but 
not all of them are equally important. Considering the 
many sources of uncertainties, the agreement between 
the predictions of conventional DE models with the 
EGRET data is remarkable: the discrepancy is only a 
factor of two at worst; this is the famous "GeV ex- 
cess." With some reasonable assumptions, such as 
CR intensity fluctuations, agreement with the data can 
be obtained in the framework of conventional astro- 
physics (21 . 

The Gamma Ray Large Area Space Telescope 
(GLAST), to be launched in 2007, will have improved 
sensitivity, angular resolution, a wider field of view, 
and uses a different technology |32| than EGRET. 
Therefore, many of the instrumental uncertainties as- 
sociated with EGRET will become irrelevant. How- 
ever, issues associated with CR propagation, gas dis- 
tribution, and other astrophysical input will remain 
waiting for the next generation of instruments. 
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